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Abstract 

The method of functional integral bosonization is extended to examine the transmis- 
sion properties of semi-infinite nanowires. In particular, it is shown that edge states 
will arise at the end point of the dimerized semi-infinite spin-chain and by combining 
the method of image and the bosonization technique, the system can be properly 
bosonized. Based on the bosonized action and a renormalization group analysis, 
it is shown that unlike scattering due to single bulk impurity in the nanowire, the 
scattering potential remains relevant even for slightly attractive potential due to the 
interaction between the edge state and its image. When the strength of potential 
goes beyond a critical strength, the tip of the semi-infinite nanowire may become 
insulating. 
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1 Introduction 



Nanowires have been considered as ideal tools for probing nano- materials [T|[2]. 
In real applications, the scanning probe microscopes (SPM) represent the typ- 
ical realization [3]. In this approach, the tip of the probe is crucial. It is there- 
fore important to understand the properties of the tip. From theoretical point 
of view, the tip can be ideally considered as the end of a one dimensional 
nanowire. The physics is thus embedded in the semi-infinite nanowires. In 
various situations, it is known that edge states may arise at the end. Never- 
theless, conventional studies of the ID wire have mostly been focused on its 
bulk properties, whereas assembled nanowires can only have finite lengths and 
must terminate at some sites (the ends, or the edges). It is therefore important 
to investigate effects that are due to the ends. 

In the absence of interactions, the presence of the edge can be handled by 
method of image or its generahzation [l|5]l6] . In this case, effects due to the 
ends are identified and reflected in the tunneling spectrum. For example, be- 
cause the end couples k and —k, the wavefunction is modulated by sin{kx) 
so that the van Hove singularity is suppressed in the local density of states 
near the end [l|5|6] . The more interesting situation occurs when the reflection 
symmetry is broken. In that case, it is shown that localized edge states may 
arise and result in peaks in the tunneling spectrum [7J. These conclusions are 
based on free electrons in which the problem is linear so that method of im- 
age is applicable. In real materials, particularly in nano-materials, interactions 
can not be neglected. In that case, the applicability of method is called into 
question. It is therefore of important to examine the applicability of method 
of image or find a general method to handle the boundary. 
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In this work, we shall show that indeed the method of image can be gener- 
alized even in the presence of interactions. The generalization results from a 
combination of method of image and the bosonization technique. By adopting 
the method of the functional bosonization, recently developed by Yurkevich[8], 
we show that the boundary conditions can be easily implemented. The mar- 
riage of the two methods allows one to investigate effects due to the edge 
states in the most general situations. To illustrate our method, we shall con- 
sider the simplest situation that supports the edge states : spinless electrons 
in dimerized semi-infinite nanowires, which would correspond to fermionized 
spin-chains. We shall first establish the existence of the edge state in the 
presence of interactions. This will be investigated by considering a dimerized 
spin chain. The method of functional bosonization is then combined with the 
method of image to calculate the two-point correlation function for dimerized 
nanowires. Finally, the bosonization of the partition function for semi-infinite 
nanowire will be derived. Based on the derived partition function, we per- 
form a renormalization group (RG) analysis and obtain the RG flow equation. 
The implication of the effect of the edge state on itinerary electrons will be 
discussed. 



2 Existence of the edge state 



We shall first illustrate that under appropriate conditions, the edge states exist 
and persist even in the presence of strong interaction. For this purpose, we start 
by considering a dimerized semi-infinite spin-chain in which the transverse spin 
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coupling changes alternatively. The Hamiltonian can be written as 

i=oo 

Hs=Y: ms:St,, + S!Sl^ + J.{St + l/2){St+, + 1/2) + hSt, (1) 

i=l 

where h is the magnetic field, Jj = ti when i is odd and Jj = ^2 when i 
is even. By applying the Jordan- Wigner transformation ?/^| = K{i)Sf and 
= S~ K{i) with -ft'(i) = exp[m YJ'f^^S^j + 1/2)], i/5 is then mapped into the 
Hamiltonian of spinless electrons with interaction between nearest neighbors 

i=oo 

H=Y. jMiJi+i + i^l^i^i) + + /^(v^J^i - 1/2). 

i=l 

Here the applied h field is equivalent to the chemical potential. In the absence 
of the Jz term, H reduces to free electron model with dimerized hopping 
amplitudes. This is the model that can describe polyacetylene[9]. In this case, 
the bulk dispersion can be solved exactly and is given by ej. = tf + t"^ + 
2tit2Cos2ka with a being the lattice constant. It is clear that the spectrum 
has a gap 2A = 2\t2 — For semi-infinite chains, previous studies show that 
when ti < t2, there is an edge state that decays exponentially from the edge 
with energy right at zero [7]. It is conceivable by continuity that the edge state 
survives when Jz is turned on as long as Jz is small in comparison to the gap. 
To check how the edge state evolves in the presence of interaction, we solve Hs 
(instead of H) by resorting to method of exact diagonalization. The energy 
of the edge state can be found by tuning h and measuring the local density 
of electrons at the 1st site. As h sweeps across the edge state, a large change 
of the local density occurs due to the fact that the edge state is localized. In 
Fig. 1, we show the numerical result of 7 sites of spin chain. The jump in the 
local density is clearly manifested. 
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The energy of the edge state can be read out by finding the magnitude of h at 
the jump. In Fig. 2, we show the energy of the edge state versus J^. Clearly, 
one sees that in presence of interactions, the edge state persists except that 
the energy is no longer right at zero. 



3 Two-point Green's functions in dimerized infinite chains 

To understand how the locahzed edge state interacts with extended states, 
we shall first consider the case without the boundary to illustrate the effect 
of dimerization on the extended states. In this case, the Hamiltonian of the 

extended states for spinless electrons is 



H = y^tiipjipi+i + h.c. + VniUi+i. 



(2) 



Here when i is odd, ti — ti, otherwise, ti — t2- In the absence of V, there are 
two bands separating by a gap. We shall set the Fermi energy £/ in the upper 
band. After the Fourier transformation, the kinetic term becomes 



^0 = E(^L(^)'CenW) 
k 



I \ ( 

Afc 







'4^even (^) 



(3) 



where = tie*'^" + t2e~'^^°'. Note that because the lattice consists of two 
sublattices: even sites and odd sites, the theory is a two-component theory 
with the obvious notations: ipodd represents -0 at odd sites and ifj^^^n represents 
■0 at even sites. To focus on the itinerant electrons, we linearize Hq in k near 
the Fermi wavevector i/c/. By replacing i{k ± kf) by dx, the action for the 
kinetic terms becomes 
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( \ 



R 



ef + ivf*dx dr — Sf 
\ 



(4) 



J 



/ 

Here Vf = a{2tit2sin2kfa — itf + it\)/ef and phase factors, e*^ = {tic^^i"- + 
t2e~'^''f°')/ef, have been absorbed into ip^^^^n V'i^en- ^ represent right 

and left movers with corresponding energy dispersions: iu = —Vfk and iuj = 
Vfk, found by diagonahzing 5*0 [10]. Note that Vf = Re[vf] is the Fermi velocity. 

In the absence of interaction, the Green's function can be found by rotating to 
a new basis so that the action is diagonalized. For example, for the left mover, 
the relation between the new basis and the old basis is (copper? ^lower)^ = 
M('?/'o ) V'i')^ with M being given by (in the k space) 



M 



V2 



exp 



■ilm[vf]k-i 
■ J 



riIm\vf]k-\ r — tirniv 

To the 1st order in k, the action for the left mover becomes 



exp[ 



-ilm[vf]k-i 
2i7 J 



-ilm[vf]k-i 



V 



(5) 



/ 



1.1 ^ 



\ 



iuj — Vfk 

iuj — 2ef + Vfk 



(6) 



Cuppers 4owerf- Obviously, 0^^^,^ and cf„^,^ correspond to bonding 



where 

and anti-bonding bands respectively. Since the lower band is filled and so we 



are left with only copper whose two-point correlation function is 
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By using Eq.(l5]), the relations between c and ip are 



jjR(x t) = -^c^ (x - ^^"-f^ t) 

Yo V-^; ' / /2 upper 2ef ' ' 

{ ■ (8) 

Ipg (x, t) — "yj copper (-^ ~l 2£j^ ' ''") 

We thus obtain the Green's functions for electrons 



/^L /^L 1 

^0^0 — '^e^e — 2^VfsinliT{(T-T')+i{x~x')/vf)//3]^ 

rL - 1 (q) 

2l3vfsinl-K((T-T')+i(x-x' jy^)/-u/)//3] 

G 



2l3vfSin[iT({T-T')+i{x-x'+ '''^^y -)/vf)/l3] 



Here except for the position shift, A = ^^^^y^, between even sites and odd 
sites, the Green's functions are the same as those for a single band. From the 
point of view of the wavefunction, if one denotes ipodd = g^^^, one can write 
^even = e^'^^'^^^^^ with 5{k) being the phase shift. Thus linearizing with respect 
io kf results in a shift of x. We then obtain 

A = dkSikf) = -a-, ~ -, (10) 

^ ^' + ti + 2tit2Cos2fc/ ^ ^ 

Therefore, the position shift can be attributed to the phase shift between even 
sites and odd sites. Clearly, when ^2 ^ ^i, A approaches —a. In this case, 
Godd->eveni^ ~ = ^odd->oddi^ — (x' + a)). This is cousistcut with the fact 
that when t2 ^ ti, electrons on two neighboring sublattice points have the 
same amplitudes. 
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We now include the effect of the interaction which takes the following form in 
the continuum approximation 



1 r 1 

-Vo = J dxdx'-{po{x),pe{x)) 



( \ 





/ 



po[x') 



pe[X 



(11) 



where p = p^ + p^ is the density operator at odd sites(o) or even sites(e) and 
Vo(x, x') = |Vo[(5(x — x' — a) + 6{x — x' + a)]. Following we first apply the 
Hubbard-Stratonovinch (HS) transformation to decouple the density operators 



e 4 



Vo 



Z 



'12) 



where Vq ^ 



V,~' 



\ ( \ 







and Z = J Dae . For fixed a 



the —iap term in combination with Eq.([5]) constitute the quadratic term. 
However, due to the matrix nature of the action, the a fields can not be 
expressed as a phase factor, e*"^, as what happens in the case when ti = t2^. 
To avoid this difficulty, we first perform the transformation of Eq.([5]). We shall 
take the left mover as a demonstration. The right mover can be handled by the 
same method. To 1st order in k, the transformed action MSkM~^ becomes 



S 



dr / dxc^^ 



Or + ivfdx — icr^ 



\ 



(13) 



where (T+ = (do + o"_)/2 and (T_ = (cTq — 0"e)/2. In the above derivation, 
we have neglected the commutator of dx and a because it is of 2nd order in 
Vo and 1//5. Furthermore, its effect is to induce a shift in the position of a; 
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equivalently, it induces a shift in the space dependence of Vo{x,x'). Since we 
shall be concerning the long distance behavior, this position shift should not 
matter. Because the anti-bonding states are 2ef below the Fermi energy, we 
can drop the anti-bonding states and keep the bonding states with the action 
being given by 

) ■ (14) 



At this stage, the effect of a can be absorbed as a phase factor by defining 



Here if we require 



{dr + iVfd^)(f)L{x, r) = a+. (16) 



the iap term is canceled out so that copper is completely free and its two-point 
correlation function of copper is exactly the same as Eq.([7]). As a result, the 
correlation function < c^pper^upper > i^ determined by < e*''^^(^''^)~*'^^(°'°) >. 
Furthermore, the relation between cr+ and (pL in Eq. f|T6|) implies that to eval- 
uate < e*'^^(^'^)~*'^^(°'°) >, one needs to find the Lagrangian for cr+. 

The Lagrangian for cr+ can be found by first noting that the transformation 
of Eq. (ITBl) induces a Jacobian 
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dr - Ef -iao Sf + ivfd^ 



InJ^ = Tr[ln- 



V 



J 



By expanding in powers of a, the above Jacobian becomes 



(17) 



n=l 



n 



where G is the 2x2 matrix of the Green's function in the absence of interac- 
tions and 



a 



( \ 

(Jn 



(19) 



Using the loop cancelation theorem [8], we retain only the n = 2 term 



(20) 



After adding contributions from left and right movers and retaining only the 
lowest orders in q and u;, we obtain 



J=J^ + J 



( \ 
1 1 



(21) 
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Combining ^ j with ^, we obtain the action for the a fields 







,2 



(22) 



Using Eqs. lITB]) and (B^ . we obtain 

< r)0i(O, 0) >= lnsin[- — {ix + vjt)] 

(JVf 

g + (yf"-*- , I n vr r r ^ / 

ln\sin[—{ix + f r)J | + iarg[sin[—{ix + f rjJJ, (23) 

where v"^ = f j(l + Vo/vrt;/) and g = Vf/v. Therefore, the two-point correlation 
function of c^pp^r is given by 



< (t rV^^ fr' r') > 

^ '-'upper V"^' ' J^upperK-^ ' ' / ^ 

1 g«a»'9[s«"[;g^(«(a;-3::')+'!'(-r--r'))]] 

(24) 



pVf {sm[n{{T - r') + i{x - x')jv)l^\)-^9'^s'^~)l'^ ' 

It is clear that the interaction changes not only the exponent but also the phase 
in correlation functions. Similar procedure applies to the right mover, the 
correlation function for the right mover turns out to be the complex conjugate 
of < c^pper^ulper >• Finally, the two-point Green's functions for electrons are 
obtained by combining Eqs. (jSj) and ([23]). The net effect of the relation Eq. 
(IE]) is to shift positions of even sites by A, given by Eq. (|TOl) . Clearly, the shift 
is not changed by the interaction and the same conclusion applies to other 
Green's functions. 

We conclude this section by noting that the bosonization of dimerized chains 
is almost the same as that for one band model except when calculating Green's 
functions, positions have to be shifted appropriately as indicating in Eq.([9]). 
Obviously, this also applies to the general case when the system contains more 
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than two sublattices. In that case, one first bosonizes the energy band cut by 
the Fermi energy. The important information one needs for calculating Green's 
functions is then contained in the transformation matrix, similar to Eq.([5]), 
which determines the position shift of each site. 



4 Generalized method of image for semi-infinite chains 

We now apply method of functional bosonization to semi-infinite chains. A 
generalized method of image will be developed in the presence of interactions. 
We first set the boundary right at site 0, thus ip{0) = 0. In the continuum 
approximation, the Hamiltonian for the energy band cut by the Fermi energy 
can be generally written a.s H = + + Vq with the interaction term being 
given by 

^ oo oo 

Vo = - y (ix y dx' p{x)Vo{x — x')p{x'), (25) 



and the kinetic terms be given by 

oo 

Hr/l = Tivf J dx^^^/^d^tpR/L- (26) 



Here the operator ip = e^^^^ipR + e~'^^^^ipL is the continuum operator for elec- 
trons in the energy band cut by the Fermi energy. In an infinite chain, the 
right-mover and the left-mover are independent. This independence, however, 
is lifted in the semi-infinite chain due to the boundary condition 

^(0)=V^^(0)+V'l(0)=0 . (27) 
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To satisfy the boundary condition, it is useful to extend the defining domain 
of ifjR and ijjL to X < For this purpose, we define 

iIjr/l{,-x,t) = -iPl/r{,x,t), (28) 

so that the boundary condition is automatically satisfied. The kinetic term 
can be then rewritten in terms of a single chiral field 

oo 

Hr + Hl = -IV f J dx^lid^ijR. (29) 

— oo 

Thus in the absence of interaction, the Green's function for the right mover is 
exactly the same as that of the infinite chain. By using Eq. (!28|) . we find that 
the full Green's function can be written as 

g{x, x'; r, r') = e'^f^'-'''^G'^{x, x'; r, r') + e-''=^("-"')G^(-x, -x'; r, r') 

_^ikfix+x')QR^^^ -a;'; r, r') + -e-''=^-("+"')G'^(-a;, x'; r, r'), (30) 

where we have denoted the Green's functions for the right/left mover by G^/^ . 
Note that Eq.([28D also implies G^{x,x') = G^{-x,-x'). Therefore, Eq.(l30l) 
can be cast in the form 

g{x, x'; r, r') = G{x, x'\ r, r') — G{x, —x'; r, r'), (31) 

where G{x,x';t,t') = e'''^'^'^^'''^G^{x,x';r,r') + e~''''f^''~''"'G^{x,x']T,T') has 
the form of the Green's function for the infinite chain. Notice that without in- 
teractions, G^/^ are the same as the bulk Green's functions for infinite chains, 
thus Eq. fl31l) demonstrates the validity of method of image. In the presence of 
the interaction, G^^^ and G appear not to be the bulk Green's functions and 
their meaning requires further exploration (see below). Nonetheless, the form 
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of Eq. (l3T]) is still valid except that one has to resort to Eq. (l30!) to calculate 
the Green's function. Obviously, all the information needed is in G^. 

To calculate G^, we first apply the HS transformation [S] on the interaction 
term for < x < oo 

^-JodrVoir) ^ J jj^^-^ J^ry;^ dx dx'aVa-\T)a+t dxac^c] ^^2) 

where Vq~^{x', x") is the interaction potential that satisfies / (ix'Vo(a;, x')Vq^{x' , x") = 
6{x — x"). Since we shall be interested in the local approximation in which 
X ~ x' dominates, by extending the HS field, o"(x, r), into — oo < x < 

a(— X, r) = cr(x, r), (33) 



and assuming that Vq^ is symmetric in the real space, the first term on the 
RHS in Eq. (!32l) can be rewritten as \ Jq dr dx dx'a{x)VQ'^a{x'). For 
the second term on the RHS in Eq. (l32l) . we use Eq.( l28l) and neglect terms 
that contain the fast-oscillating factor, e^'^^'^-f^. The action for the right mover 
then becomes 



/3 oo 

t 



Sr = / drl / dxip\{dr — i(j{x) — ivfdx)ip 



R 



-oo 

oo oo 



+ ] I dx I rfxV(x)l>o"V(x')]. (34) 
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-oo — oo 



We have thus transformed the semi-infinite problem into an infinite one except 
that the a field has to be even: cr(x) = cr(— x). Clearly, the loop cancelation 
theorem p] is still valid so that the action of a is quadratic. By making use of 
cr(g, Q) = o-{—q, Q), we obtain 
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where = f ? + Hi^. Note that the correlation function of a is given by 



(35) 



< Q') 2rr(3Vo—^-^{6{q + q') + 6{q - g')) W- (36) 

Here the term 5{q + q'), is due to the relation a{q, Vl) = cr(— fi); clearly, it is 
the crucial difference between the infinite chains and semi-infinite chains and 
reflects the effects of boundary condition. 



Similar to Eq. (fT6!) . a can be absorbed by a phase (pR with the relation 4>R{q, Q) = 
a{q, Q)/{iQ + vjq). Therefore, the correlation for the phase can be computed 
by using Eq. fl35|) . After performing appropriate contour integrations, we find 
that the analytic form for the correlation function of the phase is given by 



{(f)R{x,T)(l)R{x',T')) 



In 



sin{— — [i{x — x') — v{t — r')]} 



9 + 9~^ 



In 



TT 



sin{ — [z(x — x') — v{r — r')]} 



I arg 



TT 



sin{ — [i(x — x') — v{r — r')]} 



9 -9 



In 



TT 



sin{ — [i(x + x') - vir - r')]} 



(37) 



The Green's function of the right mover is obtained by computing e" 
we find 



■\<{<t>R{x,T)-<t>n{x' ,t')?> 



X |F(x, -X- r, r) | ^ |F(x', -x'; r', r') | ^ , 



(3J 



where F{x,x';t,t') = sin[-^{i{x — x') —v{t — t'))]. From G^, one obtains the 
full Green's function by using Eg. (1301) or Eq.( l3T]) . Finally, the Green's functions 
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of real electrons are obtained by shifting positions according to Eqs.® and 

m- 

In conclusion, the method of image, Eq. fl3ip . is valid even in the presence 
of interaction except that now the bulk Green's functions have to be re- 
interpretated as the Green's function of the infinite chains with the even a 
field. 

5 The effect of the edge state on itinerary electrons 

As demonstrated in Sec. (2), when ti < t2, a localized edge state arises. In the 
absence of interactions, its energy is right at zero with the wavefunction be 
given by ipoix) ^ (1, 0, -e, 0, e^, 0, -e^ 0, ■ ■ ■) with e = ti /to jim ?]?]: while 
in the presence of interactions, the edge state persists with the energy being 
shifted away from zero. 

Since the edge state is localized, it causes higher electron densities near the 
edge when the Fermi energy exceeds the energy of the edge state. The density- 
density interaction then induces the potential scattering on electrons near the 
Fermi energy. This is entirely similar to single bulk impurity in the nanowire, 
where if there is no interaction between the impurity and the nanowire, the 
level for the electron residing on the impurity also lies outside the energy bands 
of the nanowire. The potential scattering due to single bulk impurity was 
investigated by Kane and Fisher [TT]. In the following, instead of performing 
more detailed analysis such as calculating changes of density of states [T2], we 
shall follow Kane and Fisher and perform the stability analysis on the itinerary 
electrons qualitatively based the renormalization group analysis. The result 
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will be conceptually useful in considering the transport properties of semi- 
infinite nanowires. 

We first note that the electron operator can be represented by + V'o(x) 
with ijj{x) denote the operator for itinerary electrons near the Fermi surface. 
Substituting this representation into the interaction potential, Eq. (l25|) . a po- 
tential scattering term is induced: 

oo oo 

Vp = J '^^ J dx'Vo{x, x'){'ijjo{x'y'ipo{x'))'ip\x)'ip{x). (39) 



In the continuum approximation, we have ip{x) = e^'^f^ipR^x) + e~'^^f^ipL{x) so 
that 

oo 

% = ^f j dxv{x)[i)^i^{x)i)R{x) +^\{x)^l{x)] 



oo 

+ j dx[A(a;)4(a;)^i(x) + A(a;)*V'i(a;)^i?(x)], (40) 



where the second equality is obtained by writing in terms of ijjR and ipi with 
Vfv{x) = /o°^da;'\/o(x,x')|Vo(x')P and A(x) = /o°° dxVo(x, x')|Vo(x')| V^^^^^ By 
using Eq. fl251) and defining v{—x) = v{x) and \{—x) = A*(x), the integration 
domain can be extended to — oo. After being combined with the kinetic term, 
Eq. (l29!) . the first term on the RHS of Eq. (HO!) becomes —ivf ^^^dxip'^^ldx — 
iv{x)]ipR{x). Therefore, the first term can be absorbed into the phase of ipR 
by changing ipR ^ e^S '^^^'^^^t/^/j. The action induced by the edge state, Sedgei 
is thus determined by the second term. By writing = ipe"^^, A = |A|e*''', and 
using Eq. fl^Sl) . we find 

13 oo 

Sedge = -Jdr J dx|A(a;)|[e^''^Jj(x,r)V'ii(-x,r)e^<^«(-^'^)-*^«(^'^) 

-oo 
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where we have made use of the relation = x), which results from 

the mirror extension of the a field [see Eq.f l33p ] and the relation between a 
and the phases (pR/L [see Eq.flT6l)]. Obviously, the relevant phase is 02 (a;) = 
^{4>r{^)~(Pr{~^)): whose correlation function can be found by using Eqs. fl36l) 
and ffTB]) . We obtain 



\sin\ 


[z-z')\9 


\sin{ 




\sin{ 


>f-^'f)\ 


\sin\ 


[Z-Z')\9 



where ^ = (x, r), Zf = ^(r + ix/vf), z = ^(r + ix/v), and z is the complex 
conjugate of z. Note that terms that contain the factor Zf — z' f or z — z' are 
contributions from the image, representing the boundary effects. Following 
Ref. [12], the partition function with the edge state Z^dge =< e~^'"^!>e > can be 
expanded as 



1 /■ " 

^edge = Yl j ^ l[[dXidTi\\{Xi)\ 

n,ni+n2=n '^l-''^2- •' j 

ni ni+n2 
i j=ni+l 

(43) 



where = (— Xj, Xj). By defining (f^i = zj^i for i < rii and = for i > rii 
and applying the Wick's theorem for averaging over ipR, we find 



ni _ ni+n2 _ i 



ni ni+n2 

(44) 



On the other hand, the average of the phase field 02 can be written as 

< e2^r02{c.)-2E;i:;^i</'2fe) >^ [p(2)]7P(^^), (45) 
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where 



By using the Cauchy formula [13], one obtains det[l / sin{(f^i—Qj)] = i"''^~"'^P{zf) 
Therefore, P{zf) in Eq. fHSl) gets cancels exactly. As a result, only the velocity 
V is retained in the partition function, while terms that depend on ix ± vjt 
are absent. We obtain 



1 1 /■ " 



n,ni+n2=n 



(47) 



We now bosonize this system by showing that Zedge can be reproduced by a 
bosonic field {}. For this purpose, we first consider a free boson field, f{x), 
described by the action 

So = ^ f dxdT-{d^^Y + v{d,^ f. (48) 
Qirg J V 

The correlation function for (f is —2(7 In | sin(z — z')\. To include the effect of 
boundary, we include the image field ip{—x). The ^ is then constructed as 

^{x,r) = ^{^{x,T)-y^{-x,T)). (49) 
Hence the correlation function for is 

{'d{x,T)'d{x' ,t')) = —2g In | sin(2; — z')| — In | sin(2; — 2;')| . (50) 
It is then straightforward to verify that Z^dge can be written as 

P oo 

Zedge = {exp[— Jdr J dx\\{xW'-'^-''' + e-''-W2)])_ (51) 



-oo 
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In other words, after the bosonization, the action induced by the edge state is 



Clearly, Sedge has the same form as the action induced by impurities in a 
bulk Luttinger liquid[12]. However, the phase field d involved is different. 
As shown in Eq.( H9l) . the main difference lies in the contribution from the 
image field 0(— x, r). Following [TT], we perform a gradient expansion. First, 
since in the local approximation, V{x,x') ~ Vq5{x — x'), we have A(x) ~ 
K)|V'o(a;)pe^2'=/^ For x > 0, \Mx)\^ ~ e-2(^-«)/« with ^ = (IntsA^a and 
a being the lattice constant, X{x) peaks at x ~ a. Therefore, after taking 
A(— x) = A(x)* into consideration, in the gradient expansion, |A(x)| can be 
approximated by Ao[5(x — a) + 6{x + a)] with Aq ~ |Vo/ dx\ilJo\'^e^^^f^\. As a 
result, the action induced by the edge reduces to 



Note that rj is the phase of Vq / dxltpol'^e^'^^^^ . By setting kf ^ 1/a, rj ^ 
tan^-"^ C,/ci = tan^-'^(lnt2/^i) + 7rsi(7n(Vo). Hence when t2 = ti, we obtain 77 = 
(repulsive) or ip (attractive) and thus Sedge = 0, which is in consistent with 
the fact that the edge state disappears in this case. 

To completely describe the edge, one needs to specify the action of {}{x, r) at 
X ~ a. For this purpose, we integrate out 0(x) in Eq. (l48l) for all x except for 
X = ±a and obtain 




(52) 




(53) 
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( \ 

\ _^-2\uJn/v\a 



Aug ^ 1 - e-4|'^"/^l'=' 



_^~2\uj„/v\a 



(54) 



Here Un is the Matsubara frequency and $-^(ic<j„) = [0(a, icUn), 0(— a, -icUn)]. 
By diagonalizing 5*0, we find that the eigenmodes consist of symmetric and 
anti-symmetric combination of (j){a,iuJn) and (f){—a,iuJn)- Since only the anti- 
symmetric mode is involved in Sedg, the action for ^ is 



S^ = y^-L- 1"^"' , , \^{a,tur,)\^. (55) 



Following [TT], we perform the renormalization-group (RG) transformation. 
To leading order, the RG flow equation is 

^ = [l-^(l-e-^^«/^)]A„ (56) 

where = Aq sin rj and A is a high-frequency cutoff. Eq. fl56|) is almost the 
same as the RG flow equation for a single impurity in the bulk Luttinger 
liquid except for the correction term, e"^^"/'', due to the interaction of the 
edge state with its image. Since v = vj/g and A ^ 0{kfVf), e"^^"/"" = e"""^ 
with a being a numerical factor of order 0(1). The critical value of g below 
which A^ grows and is relevant is determined by 1 = (^£(1 — e^"^''). Clearly, the 
exact magnitude of gc depends on the microscopic details through a and is 
not universal. Nonetheless, quahtatively for Vq > TiVf{l/ g^ — 1), Vq is relevant 
and grows indefinitely under RG transformation. The tip of the nanowire thus 
may become insulating. This effect is entirely due to the interaction between 
the edge state and its image. As a result, the critical strength for the tip of 
the nanowire being insulating is slightly negative. 
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6 Conclusion and acknowledgments 

In conclusion, we have generalized the method of image to be valid in the pres- 
ence of interactions. The generalization results from a combination of method 
of image and the bosonization technique. The marriage of the two methods 
allows one to investigate effects due to the edge states in the most general situ- 
ations. Based on the derived partition function, it is demonstrated that unlike 
scattering due to single bulk impurity in the nanowire, the critical strength 
of interaction when the tip of the semi-infinite nanowire becomes insulating is 
shifted to be shghtly attractive due to the interaction between the edge state 
and its image. We gratefully acknowledge discussions with Prof. Hsiu-Hau Lin 
and the support from the National Science Council of the Republic of China. 
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Figure Captions 



Fig. 1 The local density at the 1st site versus h for zero temperature and 
l/keT = 100. Here h = 1, ^2 = 4, and = 0.1. 

Fig. 2 The shift of the edge state energy, i.e., with increasing in zero tem- 
perature. 
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